Creation of a monopole in a spinor condensate 
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We propose a method to create a monopole structure in a spin-1 spinor condensate by applying the 
basic methods used to create vortices and solitons experimentally in single-component condensates. 
We show, however, that by using a two-component structure for a monopole, we can simplify our 
proposed experimental approach and apply it also to ferromagnetic spinor condensates. We also 
discuss the observation and dynamics of such a monopole structure, and note that the dynamics of 
the two-component monopole differs from the dynamics of the three-component monopole. 
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Introduction — The experimental realization of spinor 
Bose-Einstein condensates makes it feasible to ex- 
tend the study of topological quantum objects into an 
entirely new field of physics. Ordinary single-component 
condensates have many topologically interesting proper- 
ties such as the existence of vortices But in spinor 
condensates one can also study phenomena that cannot 
exist in the single-component systems. One example is 
a monopole structure in a spinor condensate with anti- 
ferromagnetic interactions as proposed by Stoof et al. re- 
cently [Q. Other novel possibilities also exist such as 
skyrmions (|,^. With Bose-Einstein condensates a cru- 
cial aspect is not only the existence and stability of topo- 
logical structures, but also the methods for their creation 
and observation, as well as their dynamics. In this Let- 
ter we address all these aspects for a monopole structure 
in an experimentally relevant case of a multicomponent 
Bose-Einstein condensate. 

A monopole is a topological defect in a vector field. It 
is characterized by a unit vector that is radial in respect 
to some unique central point (i.e. the "hedgehog" defect). 
In spinor condensates the vector quantity could be the lo- 
cal spin of the condensed atoms , but other choices are 
also possible H,^. Monopoles have been studied theo- 
retically in two-dimensional condensates Recently 
some results for the three-dimensional case (such as den- 
sity distribution, energy and dynamics of such a defect) 
have been studied by Stoof et al. Q. They described a 
monopole created in an antiferromagetic spin-1 conden- 
sate such as ^'^Na. The monopole was characterized by a 
spinor 
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where n is the condensate density and the vector m — 
±r/r is a radial unit vector and has the spherically sym- 
metric "hedgehog" structure. Stoof et al. demonstrated 
that this particular spin texture is a unique consequence 
of the unit winding number and minimization of the gra- 
dient energy. The monopole can also be displaced from 





the center of the trap without changing their argument. 
The spinor in Eq. (|^) is non-magnetized and can be 
achieved from the single-component mean-field ground 
state, = y/n{010), with local spin-rotations. Con- 
sequently, at each position it resembles the ground state 
and thus the absence of dynamical instabilities which lead 
to domain formation |]9|,|l0| is ensured. 

Two-component monopole — One is not, however, lim- 
ited to the antiferromagnetic texture given in Eq. (|l|) 
when considering monopoles. We can alternatively map 
the vector m into an effective two-component system: 
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To ensure the stability of this texture against phase- 
separation pl| , the spin-1 condensate must have ferro- 
magnetic interactions, which makes the ^^Rb spinor con- 
densate a potential candidate. In other words, the prepa- 
ration of a monopole is not limited to the antiferroma- 
gentic ^^Na system as expected before 0], if we accept 
spinors that do not have the order-parameter space of the 
ground state. Consequently, it should be noted, that the 
texture in Eq. (||) can not be produced by local rotations 
from the ferromagnetic ground state, (p = ^/n{100). 

The static properties of the two-component monopole 
are similar to those of a three-component spinor (there 
are differences in dynamics, as we shall discuss later). 
Because it should also be easier to create experimentally, 
we focus on the two-component case although our calcu- 
lations are done for the actual three-component system. 
We note that as the spinor in Eq. (^ is not the mean- 
field ground state locally, some relaxation towards the 
true ground state is to be expected. But this requires 
spin-changing processes, that are very slow and can thus 
be ignored at timescales of interest. 

We describe the spinor condensate with a multicompo- 
nent wave function and label the components as V'm with 
the spin projection quantum number m (m = 0, ±1). 
The mean-field Gross-Pitaevskii (CP) equations are ll^] 
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-|^ilV-l) 200- 
= £^0 + Aa (2^rdV'o + l^-iPV-o + l^ilVo) (3) 
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where C = + Vtrap + A, + |^oP + l^iP), 

A, = ^(ao + 2a2), Aa = ^(aa - ao), and ap is 
the s-wave scattering length in the total hyperfine two- 
atom i^-channel. For a cylindrically symmetric trapping 
potential we use Vtrap = muj^{x^ + y^)/2 + 

Monopole structure and stability. — To understand and 
test our approach for monopole creation we study the 
monopole structure by solving Eq. (H) numerically. In 
order to create the monopole as an initial state of our 
numerical study we force the texture given in Eq. (H) (or 
Eq. (|l|)) into the order parameter (this should not be con- 
fused with the actual proposal for experimental creation 
of the monopole, to come later). The imprinted spinor is 
then propagated in imaginary time until sufficient con- 
vergence is reached. If the monopole is at the center 
of the condensate the imprint has to be done only once 
at the beginning of the iteration, otherwise the imprint 
must be repeated in the course of the iteration to prevent 
the monopole from drifting away from the intended loca- 
tion, to a location with lower energy. In Fig. |^ we show 
the typical density distribution of the spin-1 monopole 
located at the center of a trap. 

By looking at the individual condensate components 
we gain relevant insight into the structure of the 
monopole. Here tp-i has a vortex at 2; = with a core 
size that is a function of z. On the other hand, ipo goes 
through a tt phase shift as we move from positive to neg- 
ative z values; consequently, this component relates to 
a soliton in a single-component condensate. The m = 
component atoms fill the vortex line everywhere else ex- 
cept at the origin, where the density of m = compo- 
nent also vanishes. Therefore the intersection of the vor- 
tex line with the soliton plane gives rise to a monopole 
core. Thus, if we can experimentally create a vortex and 
a soliton in a two-component system, we can obtain a 
monopole. 

If we displace the monopole from the center of the 
spherically symmetric trap, the energy of system de- 
creases as function of the displacement (in qualitative 
agreement with the results in Ref. Q). This indicates 
that in a dissipative environment the monopole is ex- 
pected to drift to the edge of the condensate and van- 
ish. The estimate for the timescale of this process goes 
beyond the model used in this paper. But in single- 
component condensates the vortex lifetimes can be sev- 
eral seconds jl^ and this timescale is also, presumably, 
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FIG. 1. The total density {I = \/h/muj) of the spin- 1 
condensate with 5- 10* *^Rb atoms when y = z = 0.05 (origin 
is not present in our discretization scheme). The monopole 
is at the center of a spherically symmetric trap with a trap 
frequency uj = ujr = i^z = (27r) 50 Hz. 

indicative of the monopole lifetime. On the other hand, 
the dissipative processes might be masked by the topo- 
logical instability of the monopole. The unit winding 
number of the spin texture is not sufficient to protect the 
topological stability of the system and instead of drifting 
smoothly to the condensate edge, the monopole might 
decay into other excitations. 

Creation of monopoles. — The separate look into each 
spin component of the monopole structure suggests a pos- 
sible way to create it. For example, we can prepare a 
spinor condensate with 2/3 population at the m = 1 state 
and the rest at m = state, e.g. with an rf-pulse [ p^ . 
A "blueprint" of the monopole is achieved by creating a 
vortex into the m = 1 component and a soliton (with 
TT phase discontinuity) into the m — Q component, both 
at the trap center. The vortex line should lie along the 
phase discontinuity in the m — Q wavefunction. In Fig. ^ 
we demonstrate the time-evolution of such a mixture in a 
cigar shaped trap in real time, when the soliton and vor- 
tex were created using phase- imprint method . With 
other excitations abound, it is clear that we nevertheless 
have a monopole inside the condensate. 

In our numerical studies we have used a certain amount 
of smoothing to reduce the amount of noise that would 
be created if a phase-imprint is too abrupt. Smoothing 
for a vortex was done by assuming that not only do we 
have a phase-mask, but also a narrow (on the order of the 
coherence length) laser beam that bores a hole through 
the m = 1 component along the vortex line. For a soli- 
ton the TT phase jump was done at the distance of the 
order of the coherence length. Without such smoothing 
our numerical approach becomes unstable. It is not clear 
how much is this smoothing required in actual experi- 
ments, although some amount would be always present. 
To obtain this combined phase-imprint experimentally is 
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probably complicated, but at least the two main ingre- 
dients, namely experimental creation of vortices [l^ and 
solitons has already been achieved. 
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FIG. 2. The time evolution of the total density of a con- 
densate with 5 • 10^ ^'^Rb atoms in a cylindrically symmetric 
trap with frequencies u) = ujr = (27r) 250 Hz and uj^ = (27r) 50 
Hz. Initially there is a vortex at the m = —1 component 
(perpendicular line) and a soliton at the m = component 
(horizontal line). The vortex and the soliton are imprinted 
at t = and the figures show the cut of the total density in 
y — 0.09 plane. a;-axis is along the horizontal direction. 

Observation of monopoles. — The monopole core has 
roughly the size of the healing length, and it is inside the 
condensate. Thus its direct observation is difficult. But 
the same methods used to observe vortices and 
vortex rings [ p^ can be applied to observe monopoles as 
well. One should first let the condensate expand and then 
image the 3D structure of the different m-states using two 
orthogonal probe beams. As for separating the different 
m-states, one can use an appropriate Stern-Gerlach ap- 
paratus 1^. 

We have studied the behavior of a freely expanding 
monopole using the time-dependent generalization of the 
monopole GP-equation in Ref. [Q. This approximation 
assumes equal scattering lengths, but at the timescales of 
interest the role of differing scattering lengths is in fact 
negligible. In the limit that 02 = ao all the components 
feel the same spherically symmetric potential (external 
potential and the mean field terms) and are not sensitive 
to the phase of the other components. Therefore a single- 
component GP-equation is sufficient for modelling the 
expanding monopole. 

In Fig. ^ we show an example of the time-evolution 
of the ratio of the monopole core size to the system ra- 
dius once the trapping potential is turned off. Expansion 
is qualitatively similar to vortex expansion in a scalar 



condensate. At small times the monopole core size, ^, 
adjusts (almost) instantaneously to the local density |l^] 
and one expects the size of the core to scale as the healing 
length, ^0 = l/Vdvran. If we model the wave function as 

^(r) = Acxp exp (^^^ , (4) 

where A is the normalization factor, the monopole grows 
faster than the expanding condensate, or more precisely 

This happens as long as the characteristic time for ad- 
justment of the core size, Tad ^ fi/nXs, is much less than 
the expansion time t^x ~ R/cg, where Cs is the sound 
velocity. The parameters in Fig. ^ imply that in this 
regime the condensate can expand by an order of mag- 
nitude, and in the end of this regime the monopole size 
compared to the condensate size has greatly increased. 
At later times the atoms in the cloud will evolve as free 
particles and io/R will settle to some constant value. We 
also compared the results obtained with a single compo- 
nent GP-equation (at small times) against the solution 
of the multi-component GP-equations (|^) and found that 
the two approaches give essentially the same results. 

0.2| > > 1 
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FIG. 3. Time evolution of the ratio of the monopole core 
size Rm to the condensate size Rc. These sizes were deter- 
mined from locations where the density was one half of the 
maximum density. Initially we assumed 5 • 10* condensed ru- 
bidium atoms in a spherically symmetric trap with a trap 
frequency uj — {2-k) 50 Hz. 

Monopole dynamics — The dynamics of a monopole is 
quite interesting. As expected, the monopole at the ori- 
gin is stable and stationary. A displaced monopole, on 
the other hand, behaves differently |^^. The monopole 
precesses around the trap center inside the condensate, 
just like a displaced vortex line does. It returns to its 
initial location after T = 2n/Vlp. For a vortex close to 
the center the precession frequency 17 p coincides with the 
frequency of the anomalous mode w^. For a disc-shape 
trap an analytic result is available [pl| and is given by 
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where R is the radius of the system and ^ is the healing 
length. Even though the trap geometry in our example is 
nowhere near the disk-shape we expect that Eq. (||) gives 
a reasonable order of magnitude estimate. Especially so 
since the m — —1 atoms are "squeezed" between lobes of 
m = atoms, thus making the disc-shaped approxima- 
tion rather justified. 

As a test case we take a spherically symmetric trap 
with trap frequency cu — (27r) 50 Hz and a *^Rb con- 
densate with 5 • 10'' atoms. Setting R to the Thomas- 
Fermi radius of the system and calculating ^ from the 
Thomas-Fermi result at the trap center we get an esti- 
mate ujT « 35. This value is fairly close to the value 
actually seen in our 3D simulation of the monopole dy- 
namics. If the monopole was displaced by one-fifth of the 
Thomas- Fermi radius our numerical result is loT « 38. 
The above estimate is surprisingly accurate. In partic- 
ular because the precessing vortex should feel the mean 
field of the other component. 

The dynamical behavior of the three-component 
monopole of Eq. (|^) is different from above. In this case 
one has a vortex at the m = —1 state and an anti- vortex 
at the m — 1 state. As a displaced vortex and an anti- 
vortex precess in opposite directions, the monopole core 
will vanish only to reappear at the opposite side as soon 
as the vortices have precessed that far. This recurrence 
is almost perfect [Q. Partial revival of a monopole has 
also been predicted in case of a 2D monopole Q. As 
the vortices precess in opposite directions the order pa- 
rameter becomes magnetized and can no longer be rep- 
resented in the form given by Eq. (|l|). Therefore, the or- 
der parameter is no longer in the order parameter space 
of the ground state. Obviously, our proposed approach 
to create monopoles experimentally applies to the three- 
component case as well. But then one needs to create a 
vortex/antivortex in the m = ±1 state, respectively, in 
addition to the soliton in the m = state. 

In Ref. jj] the dynamics of the monopole were due 
to the two different scattering lengths. In our inhomo- 
geneous spinor condensate the dynamics are not inti- 
mately connected with differing scattering lengths. Set- 
ting a2 = flo does not change our results qualitatively 
and even the quantitative changes are small. At longer 
times some population dynamics can be observed, but 
the dynamics of the total density is almost unaffected. 
Therefore it seems clear that in an inhomogeneous spinor- 
condensate the dynamical behavior of the monopole goes 
beyond the model suggested by Stoof et al. in Ref. Q. 

To summarise, we have proposed and demonstrated 
numerically a method to create monopoles in three- 
dimensional Bose-Einstein condensates, and shown that 
monopole creation is not limited to the antiferromagnetic 
spinor condensates. In addition we have studied the de- 



tection by expansion of such monopoles, and also the 
dynamics of displaced monopoles in a trap. For the cre- 
ation of two-dimensional monopoles, an off-resonant Ra- 
man beam has been suggested 0, but so far there has 
not been any suggestions for creating a three-dimensional 
monopole in a realistic experiment. Also, by replacing 
a single vortex at one component with a lattice of vor- 
tices in our approach could lead to creation of multiple 
monopoles, and allow one to study interactions between 
monopoles. 
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